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Abstract 

J  A  set  of  linear  programming  test  problems  is  analyzed  with  MINOS,  Version  5.1.  The 
problems  have  been  run  with  different  options  for  scaling  and  partial  pricing  to  illustrate 
the  effects  of  these  options  on  the  performance  of  the  simplex  method.  The  results  indicate 
that  the  different  options  can  significantly  improve  or  degrade  the  performance  of  the 
simplex  method,  and  that  these  options  must  be  chosen  wisely. 

For  each  problem,  a  picture  of  the  nonzero  structure  of  the  matrix  A  is  also  presented 
so  that  the  problems  can  be  classified  according  to  structure,  fp  / _ 
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2.  Testing  Methodology  and  Results  1 

1.  Introduction 

The  recent  development  of  new  algorithms  for  linear  programming  has  sparked  some  serious 
computational  testing  of  these  algorithms,  with  comparisons  often  being  made  against  the 
simplex  method  (Dantzig,  1963),  as  implemented  in  some  version  of  MINOS  (Murtagh  and 
Saunders,  1980, 1983,  1987).  For  such  tests,  MINOS  is  usually  run  with  certain  options  set 
at  their  default  values.  Changing  the  settings  of  these  options  affects  the  performance  of 
the  simplex  method,  and  hence  can  affect  the  conclusions  one  makes  when  comparing  the 
performance  of  MINOS  to  another  linear  programming  code.  The  purpose  of  this  paper  is 
to  analyze  how  much  some  of  these  options  can  change  the  speed  of  MINOS  when  solving 
linear  programs. 

Gay  (1985)  has  made  a  set  of  linear  programming  test  problems  publicly  available  via 
netlib  (Dongarra  and  Grosse,  1987).  The  form  of  each  problem  is 

•  T 

mm  c x 
subject  to  Ax  o  6, 

/  <  x  <  u, 

where  the  symbol  o  indicates  that  each  row  af  is  related  to  6,-  by  a  relation  of  the  form 
ajx  <  bi ,  ajx  >  bi ,  or  ajx  =  At  the  time  of  this  study,  53  problems  were  in  the  library. 
Table  1.1  lists  the  problems,  ordered  by  the  number  of  nonzeros.  Each  problem  was  given 
a  number  by  Gay,  which  is  reflected  in  the  last  column.  For  ease  of  presentation  of  some 
graphs  later  in  this  report,  the  problems  have  been  renumbered  from  1  to  53,  as  listed  in 
the  first  column.  The  optimal  objective  values,  as  reported  by  Gay,  are  also  listed.  Later, 
it  will  be  interesting  to  note  the  slight  differences  between  these  values  and  the  values 
obtained  by  MINOS  when  the  scaling  and  partial  pricing  options  are  varied.  Table  1.2 
gives  some  statistics  on  each  problem  indicating  the  number  of  different  types  of  rows  and 
columns.  A  normal  row  or  column  has  either  a  lower  or  upper  bound.  A  bounded  row  or 
column  has  distinct  lower  and  upper  bounds.  A  fixed  row  or  column  has  equal  lower  and 
upper  bounds,  while  a  free  column  has  no  lower  or  upper  bound. 

In  Section  2,  the  testing  methodology  is  explained  and  the  results  are  presented  for  six 
different  combinations  of  MINOS  options.  Section  3  contains  an  analysis  and  discussion  of 
these  results.  Section  4  presents  pictures  of  the  nonzero  structure  of  the  constraint  matrix 
for  each  of  the  problems. 

2.  Testing  Methodology  and  Results 

The  computational  testing  in  this  report  was  performed  on  a  DEC  VAXstation  II  with 
13  megabytes  of  main  memory.  The  operating  system  was  VMS,  version  4.5,  and  the 
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VMS  FORTRAN  compiler,  version  4.6,  was  used  with  the  default  options,  including  code 
optimization.  Each  run  was  made  as  a  batch  job,  and  the  account  used  to  run  the  tests 
was  given  access  to  all  of  the  main  memory.  Hence,  the  memory  used  for  each  test  was 
never  paged  onto  the  disk,  and  the  amount  of  paging  was  kept  to  a  minimum.  The  VMS 
timing  mechanism  used  is  slightly  affected  by  disk  input  and  output  (such  as  printing  the 
iteration  log),  and  hence  all  timing  results  contain  a  slight  but  negligible  error. 

For  each  test  run,  the  total  time  of  the  run  and  the  solution  time  were  recorded.  The 
solution  time  was  measured  by  timing  the  MINOS  subroutine  M5SQLV.  MINOS,  Version  5.1, 
(June  1987)  was  used.  In  the  course  of  running  the  experiments,  a  degenerate  cycle  of 
6  was  discovered  in  one  of  the  runs  on  the  problem  PILOT J A,  with  this  cycle  occurring 
in  Phase  2  after  43000  total  iterations.  MINOS  was  modified  to  avoid  such  cycles  and 
all  of  the  problems  were  run  again.  The  main  changes  required  were  in  the  step-length 
procedure  M5CHZR.  These  changes  will  be  documented  elsewhere  and  made  available  in  a 
later  version. 

MINOS  requires  a  specifications  file  for  each  run.  The  file  used  in  all  cases  was  of  the 
form: 


BEGII  SPECS  file  for  timing  LP  problems 

ROUS 

2400 

C0LUMVS 

10000 

ELEMEITS 

45000 

SUMMARY  FILE 

9 

MPS  FILE 

10 

ITERATIOVS 

100000 

PRIFT  LEVEL 

0 

SCALE  OPTIOV 

a 

PARTIAL  PRICE 

P 

EID 


The  options  changed  were  the  scaling  option  a  and  the  partial  pricing  option  p.  The  default 
values  are  a  =  0  and  p  =  1.  The  six  combinations  used  were  a  =  0  with  p  =  1  and  10,  and 
s  =  2  with  p  =  1,  10,  20,  and  40.  The  effects  of  partial  pricing  can  be  analyzed  with  the 
runs  corresponding  to  a  =  2,  while  the  effects  of  scaling  can  be  analyzed  when  p  =  1  and 
p  =  10.  It  should  be  noted  that  most  researchers  report  the  results  for  a  =  0  and  p  =  1, 
the  usual  default  values.  (In  some  cases,  the  default  value  for  p  may  not  be  1.) 

By  default,  MINOS  chooses  a  triangular  initial  basis  from  the  matrix  [A  /],  as  opposed 
to  an  initial  basis  consisting  of  only  the  slack  variables.  Certainly,  a  different  scheme  used 
to  choose  the  initial  basis  would  produce  different  computational  results. 
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Tables  2.1-2.6  contain  the  computational  results  for  the  53  problems  using  the  options 
above.  The  column  for  “Phase  1  Its.”  is  the  number  of  iterations  needed  to  find  the  first 
feasible  point.  The  column  for  “Phase  2  Its.”  is  the  number  of  iterations  needed  after  the 
first  feasible  point  has  been  found.  The  column  for  “Total  Its.”  is  the  sum  of  the  preceding 
two  columns.  The  column  for  “Blocked  Its.”  is  the  number  of  simplex  iterations  for  which 
movement  was  blocked,  i.e.,  the  number  of  iterations  for  which  the  MINOS  routine  M5CHZR 
returned  the  variable  MOVE  equal  to  .  FALSE .  This  happens  when  a  small  step  is  taken  and 
provides  a  measure  of  the  degree  of  degeneracy  of  each  problem.  The  column  for  “Total 
Time”  is  the  number  of  CPU  seconds  that  MINOS  used  from  beginning  to  end  (including 
reading  in  the  problem  data  and  writing  out  the  solution).  The  column  for  “Solution 
Time”  is  the  number  of  CPU  seconds  used  by  the  call  to  the  routine  M5S0LV. 


3.  Analysis  of  Results 

Since  linear  programs  are  solved  in  finite-precision  arithmetic,  the  “optimal”  objective 
value  may  change  depending  on  the  chosen  options.  Using  the  optimal  value  for  SCALE 
OPTIOV  2  and  PARTIAL  PRICE  1  for  comparison,  Table  3.1  indicates  the  relative  error  of 
the  optimal  objective  value  reported.  Each  column  represents  one  scaling  and  pricing 
combination,  and  only  nonzero  relative  errors  are  given.  The  first  column  contains  the 
optimal  objective  value  used  for  comparison. 

To  compare  the  efficiencies  of  the  different  options,  performance  factors  are  used.  Let 
t»,p  be  the  solution  time  for  a  specific  problem  when  the  MINOS  options  SCALE  OPTIOM  s 
and  PARTIAL  PRICE  p  are  used.  Then  the  performance  factor 


ft, Pi 
J»\,pi 


measures  the  efficiency  of  the  option  combination  (S21P2)  relative  to  the  option  combination 
(si,Pi).  For  example,  =0.5  indicates  that  when  full  pricing  is  used,  scaling  caused 
the  problem  to  be  solved  in  half  as  much  time  as  when  scaling  was  not  used. 

Performance  factors  can  be  computed  for  each  problem  to  compare  any  scaling/pricing 
pairs.  Figure  3.2  demonstrates  the  range  of  the  performance  factors  /J’f  for  the  different 
combinations  used.  The  horizontal  axis  is  ordered  by  problem  number,  and  hence,  by  the 
number  of  nonzeros.  Note  that  a  dot  representing  the  value  of  /J’j  (for  the  default  com¬ 
bination  of  options)  is  implicitly  plotted  at  the  value  1.0  for  each  problem.  All  dots  below 
the  value  1.0  indicate  a  combination  that  performed  better  than  the  default  combination 
for  MINOS,  while  dots  above  the  value  1.0  indicate  combinations  that  performed  worse 
than  the  default  combination.  It  is  interesting  to  note  that  the  factors  range  from  0.1 
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to  3.8,  indicating  that  changing  the  options  improved  the  performance  of  MINOS  by  as 
much  as  a  factor  of  10  and  worsened  the  performance  by  as  much  as  a  factor  of  3.8.  Each 
vertical  line  represents  one  problem,  and  the  ordering  of  the  factors  of  the  six  different 
options  varies  among  the  problems,  indicating  the  difficulty  of  choosing  a  "best”  option. 
Except  for  problems  31  and  46  (SCSD6  and  SCSD8),  the  trend  seems  to  indicate  that  par¬ 
tial  pricing  is  more  effective  as  the  number  of  nonzeros  increases.  However,  for  those  two 
problems,  the  performance  deteriorated  when  scaling  was  used. 

3.1.  Effects  of  Scaling.  In  order  to  compare  the  effects  of  scaling,  the  values  and 
/o’io  have  been  plotted  in  Figure  3.3.  The  plot  order  is  determined  by  sorting  on  /o’/,  to 
make  any  trend  more  obvious.  This  graph  indicates  that  the  effects  of  scaling  are  usually 
independent  of  the  partial  pricing  option  chosen.  (For  problem  49,  PILOTJA,  the  widely 
separated  values  are  due  to  the  poor  performance  when  scaling  was  not  used.)  Also,  71  out 
of  106  of  the  tests  exhibited  a  value  of  less  than  1.0,  indicating  that  scaling  improved 
the  performance  of  MINOS  67%  of  the  time.  50%  of  the  tests  exhibited  a  value  less  than 
0.9,  while  84%  of  the  tests  exhibited  a  value  less  than  1.1.  This  means  that  for  half  the 
problems,  scaling  gave  a  significant  improvement,  while  for  most  problems,  scaling  was 
either  better  than  or  comparable  to  no  scaling. 

3.2.  Effects  of  Partial  Pricing.  The  effects  of  partial  pricing  are  analyzed  in  Figure  3.4 
by  plotting  the  values  for  //’f  for  p  =  10,  20,  and  40.  The  problems  are  sorted  by  //'j10 
in  order  to  study  whether  partial  pricing  improves  performance  uniformly.  Here,  we  can 
see  that  partial  pricing  does  reduce  the  CPU  time  for  most  of  the  problems.  Of  course,  if 
the  goal  were  to  choose  the  value  of  p  that  minimizes  CPU  time,  the  results  indicate  that 
determining  that  value  may  require  extensive  experimentation,  even  for  a  single  problem. 

3.3.  Iteration  Counts.  The  effects  of  the  different  options  on  iteration  counts  are 
analyzed  in  Figure  3.5  by  plotting  the  values  of  a  =  I/m ,  where  I  is  the  number  of 
iterations  and  m  is  the  number  of  rows.  The  problems  are  ordered  by  the  number  of  rows 
to  reflect  any  trend  in  a  due  to  row  size.  A  logarithmic  scale  was  used  to  demonstrate 
the  range  of  values  of  a,  which  varies  from  0.2  to  50.  Thus  there  is  a  large  variance  in  a, 
but  there  is  no  apparent  upward  or  downward  trend.  In  35%  of  the  tests,  partial  pricing 
reduced  the  number  of  iterations  as  compared  to  full  pricing.  Intuitively,  one  would  expect 
that  partial  pricing  would  increase  the  number  of  iterations  used  by  the  simplex  method, 
since  the  "best”  (in  terms  of  minimal  reduced  cost)  column  is  not  chosen  on  each  iteration. 
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4.  Pictures  of  Problems 

In  order  to  develop  a  further  understanding  of  the  problems  in  this  set,  the  nonzero  struc¬ 
ture  of  each  constraint  matrix  is  presented  here.  Free  rows  and  columns  for  fixed  variables 
were  eliminated.  The  pictures  were  produced  on  an  Apple  Macintosh  Plus,  using  software 
developed  by  the  author,  and  then  printed  on  an  Apple  LaserWriter,  which  has  a  resolution 
of  300  dots  per  inch.  In  order  to  fit  the  problems  on  a  page,  it  was  necessary  to  print  the 
problems  at  various  magnifications.  A  magnification  of  1.0  represents  75  rows/columns 
per  inch,  so  that  a  4x4  box  of  pixels  on  the  laser  printer  represents  one  nonzero.  Hence, 
a  picture  at  a  magnification  of  0.25  has  300  rows/columns  per  inch,  and  every  nonzero 
corresponds  to  one  LaserWriter  pixel.  The  various  magnifications  and  the  correspondence 
between  nonzero  elements  and  LaserWriter  pixels  are  summarized  in  Table  4.1.  The  or¬ 
dering  of  the  rows  and  columns  is  determined  by  their  ordering  in  the  MPS  decks  that  are 
created  by  the  program  EMPS,  available  from  netlib. 

Magnification  Rows/Columns  Nonzeros  LaserWriter 

Per  Inch _ Pixels 

10.0 _ X5 _ 1 _ 40 

4.0 _ 18J5 _ 1 _ 16 

2.0 _ 3X5 _ 1 _ 8 

1.5 _ 50 _ 1 _ 6 

1.0 _ 75 _ 1 _ 4 

0.5 _ 150 _ 1 _ 2 

0.25 _ 300 _ 1 _ 1 

0.125 _ 600 _ 2 _ 1 

0.0625 _ 1200 _ 4 _ 1 

0.05  1500  I  5  I  1 

Table  4.1  Magnification  Correspondences 

It  is  interesting  to  note  that  a  large  number  of  problems  in  this  test  set  have  stair¬ 
case  structure.  Various  subsets  of  the  netlib  problems  have  been  used  to  compare  the 
simplex  method  with  interior-point  algorithms  (Karmarkar,  1984).  The  more  favorable 
results  reported  for  the  interior-point  approach  tend  to  be  associated  with  strong  staircase 
structure  (see,  e.g.,  Adler,  Resende  and  Veiga,  1986;  Monma  and  Morton,  1987).  This  is 
fortuitous,  since  staircase  problems  have  long  been  viewed  as  unusually  difficult  for  the 
simplex  method  (Fourer,  1982).  These  problems  tend  to  require  many  simplex  iterations 
to  solve  and  to  have  rather  dense  basis  factorizations.  It  remains  to  be  seen  whether  most 
problems  of  interest  in  the  “real  world”  display  staircase  structure. 
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5.  Conclusions 

From  these  results,  it  is  apparent  that  scaling  and  partial  pricing  can  improve  the  perfor¬ 
mance  of  the  simplex  method  most  of  the  time.  However,  when  comparing  the  performance 
of  the  simplex  method  to  some  other  algorithm,  care  must  be  taken  when  choosing  these 
options  in  order  not  to  degrade  the  performance  of  the  simplex  method. 

The  pictures  of  the  problems  reveal  similarities  between  groups  of  problems,  and 
indicate  that  the  test  set  in  netlib  may  have  some  bias  which  affects  these  and  other 
computational  results. 
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CAPRI 

SCAGR25 


SCTAP1 

BRANDY 

ISRAEL 


ETAMACRO 

SCFXM1 

GROW7 


BANDM 

E226 

STANDATA 


SCSD1 

GFRDPNC 

BEACONFD 


STAR 

SCRS8 

SEBA 


SHELL 

PILOT4 

SCFXM2 


SCSD6 

GROW15 

SHIP04S 


FFFFF800 

GANGES 

SCFXM3 


SCTAP2 

GROW22 

SHIP04L 


PILOTWE 

SIERRA 

SHIP08S 


SCTAP3 

SHIP12S 

25FV47 


SCSD8 

NESM 

CZPROB 


PILOTJA 

SHIP08L 

SHIP12L 


80BAU3B 

PILOT 


Rows 


28 

57 


CoIbmos  Nooteroes 


14706 

17085 

21597 


29063 

43220 


WTTA 


■4.64753 1428600CE  . 
225494963 16000E+05 
-5220206121 1707E+01 


-2.33 1 3892547843E- 
-4.157322407400CE 
-2.6661600000000E 


1 2983 1462461 36E  +05 
-7.65893 18579000E+04 
1.3730803942085E+03 


1 578124822738 1E+03 
2.6900129138000E+03 
-1.4753433060769E+O7 


1.4122500000000E+03 
15185098965000E+03 
-8. 9664482 186000E+05 


-7557 152278500QE+02 
1 .8416759028 349E+04 
-4.778781 1814712E+07 


-15862801845000E+02I 

-1.8751929066000E+01 

12576995000000E+03 


8.6666666743334E+00 
6.902235999 5488E+06 
3.3592485807 000E +04 


•251266951 19000E 
9.0429695380079E+02 
1.571 1600000000E +04 


1 208825 3460000E+09 
-2581016225338 1E+03 
3.6660261 564999E +04 


5.0500000078262E+01 
-1.0687094129358E+08 
1 .7987 147004454E+06 


5556799608 5000E+05 
-5247358 1036000E- 12 
5.490125454975 1E+04 


1 .724807 142857 1E+03 
-1.608 3433648256E +08 
1.7933245379704E+06 


-2.7190896937050E 
1.5394362183632E+07 
1.9200982105346E 


1.4240000000000E+03I 
1.4892361 344061E+06J 
5.50I8458883000E+03 


9.0499999992546E+02 

1.4076121092176E+O7 

2.1851966989000E+06 


-6.1 100885456374E+03| 
1.90905521 13891E+06 
1 .4701 879 193293E 


9.87224 19288205E+05 
-5.5742017351000E+02 


Prob.  #  (aetlib 


Tabic  1.1  Problen  Statistics  I 
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Problem 

Rows 

Name 

Fixed 

Free 

Fixed 

AFRO 

19 

8 

0 

32 

0 

0 

0 

ADLOTLE 

41 

15 

0 

97 

0 

0 

0 

SC205 

114 

91 

0 

203 

0 

0 

0 

SCA0R7 

45 

84 

0 

140 

0 

0 

0 

SHARE2B 

83 

13 

0 

79 

0 

0 

0 

RECIPE 

24 

67 

0 

85 

0 

26 

69 

VTPBASE 

143 

55 

0 

87 

1 

18 

97 

SHARE  IB 

28 

89 

0 

225 

0 

0 

0 

BORE3D 

19 

214 

0 

302 

0 

1 

12 

SCORPION 

108 

280 

0 

358 

0 

0 

0 

CAPRI 

129 

142 

0 

192 

14 

16 

131 

SCAGR2S 

171 

300 

0 

500 

0 

0 

0 

SCTAP1  j 

180 

120 

0 

480 

0 

0 

0 

BRANDY 

54 

166 

0 

249 

0 

0 

0 

ISRAEL 

174 

0 

0 

142 

0 

0 

0 

ETAMACRO 

128 

272 

0 

426 

0 

82 

180 

SCFXM1 

143 

187 

0 

457 

0 

0 

0 

GROW7 

0 

140 

0 

21 

0 

0 

280 

BANDM 

0 

305 

0 

472 

0 

0 

0 

E226 

190 

33 

0 

282 

0 

0 

0 

STANDATA 

199 

160 

0 

955 

0 

16 

104 

SCSD1 

0 

77 

0 

760 

0 

0 

0 

GFRDPNC 

68 

548 

0 

834 

0 

0 

258 

BEACONFD 

33 

140 

0 

262 

0 

0 

0 

STAR 

147 

209 

0 

373 

6 

82 

6 

SCRS8 

106 

384 

0 

1169 

0 

0 

0 

SEBA 

1 

507 

7 

521 

0 

0 

507 

SHELL 

2 

534 

0 

1399 

0 

250 

126 

PILOT4 

123 

287 

0 

635 

88 

30 

247 

SCFXM2 

286 

374 

0 

914 

0 

0 

0 

SCSD6 

0 

147 

0 

1350 

0 

0 

0 

GROW15 

0 

300 

0 

45 

0 

0 

600 

SHIP04S 

48 

354 

0 

1458 

0 

0 

0 

FFFFF800 

174 

350 

0 

854 

0 

0 

0 

GANGES 

25 

1284 

0 

1277 

0 

0 

404 

SCFXM3 

429 

561 

0 

1371 

0 

0 

0 

SCTAP2 

620 

470 

0 

1880 

0 

0 

0 

GROW22 

0 

440 

0 

66 

0 

0 

880 

SHIP04L 

48 

354 

0 

2118 

0 

0 

0 

PILOTWE 

139 

583 

0 

2335 

80 

78 

296 

SIERRA 

699 

528 

0 

0 

0 

20 

2016 

SHIP08S 

80 

698 

0 

2387 

0 

0 

0 

SCTAP3 

|  860 

620 

0 

2480 

0 

0 

0 

SHIP12S 

106 

1045 

0 

2763 

0 

0 

0 

25FV47 

305 

516 

0 

1571 

0 

0 

0 

SCSD8 

0 

397 

0 

2750 

0 

0 

0 

NESM 

94 

480 

88 

1009 

0 

175 

1739 

CZPROB 

39 

890 

0 

3294 

0 

229 

0 

PILOTJA 

279 

661 

0 

1250 

88 

311 

339 

SHIP08L 

80 

698 

0 

4283 

0 

0 

0 

SHIP12L 

106 

1045 

0 

5427 

0 

0 

0 

80BAU3B 

2262 

0 

0 

6244 

0 

498 

3057 

PILOT 

1208 

233 

0 

2320 

0 

203 

1129 

Table  1.2  Problem  Statistics  II 


I 

I 

l 

k  ^ 

5 


mm 


Problem  Nam 


AFIRO 

ADUTTLE 

SC205 


SCAGR7 

SHARE2B 

RECIPE 


VTPBASE 
SHARE  IB 
BORE3D 


SCORPION 

CAPRI 

SCAGR25 


SCTAP1 

BRANDY 

ISRAEL 


ETAMACRO 

SCFXMI 

GROW7 


BANDM 

E226 

STANDATA 


SCSDI 

GFRDPNC 

BEACONFD 


STAIR 

SCRS8 

SEBA 


SHELL 

PILOT4 

SCFXM2 


SCSD6 

GROW15 

SHIP04S 


FFFFF800 

GANGES 

SCFXM3 


SCTAP2 

GROW22 

SHIP04L 


PILOTWE 

SIERRA 

SHIP08S 


SCTAP3 

SHIP12S 

25FV47 


SCSD8 

NESM 

CZPROB 


PILOTJA 

SHIP08L 

SHIP12L 


80BAU3B 

PILOT 


Phase  1  Its.  IPhase  2  It 


3 

8 


Table  2.1  Results  for  SCALE  OPTION  0  and  PARTIAL  PRICE  1 


ms&sim 


r-A  .v  ,  »  .N 


V-L-o' 
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Problem  Nai 


ADUTTLE 

SC205 


SCAGR7 


IKSTEvkW; 


RECIPE 


VTPBASE 
SHARE  IB 
BORE3D 


SCTAPI 

BRANDY 

ISRAEL 


ETAMACRO 

SCFXM1 

GROW7 


BANDM 


STANDATA 


SCSD1 

GFRDPNC 

BEACONFD 


STAIR 

SCRS8 

SEBA 


SHELL 

PILOT4 

SCFXM2 


SCSD6 

GROW1S 

SHIP04S 


FFFFF800 

GANGES 

SCFXM3 


SCTAP2 

GROW22 

SHIP04L 


PILOTWE 

SIERRA 

SHIP08S 


SCTAP3 

SHIP12S 

25FV47 


SCSD8 

NESM 

CZPROB 


PILOTJA 

SHIP08L 

SHIP12L 


80BAU3B 

PILOT 


Phase  1  Its.  IPhase  2  ItsJ  Total  Its 


3 
17 


■EEHEMEB 


Total  Time!  Solution  Time 


9114.7 

93384.7 


Table  12  Results  for  SCALE  OPTION  0  and  PARTIAL  PRICE  10 


AFIRO 

ADUTTLE 

SC205 


SCAGR7 


IEIQ7UMT 


RECIPE 


VTPBASE 
SHARE  IB 
BORE3D 


SCORPION 

CAPRI 

SCAGR25 


SCTAP1 

BRANDY 

ISRAEL 


ETAMACRO 

SCFXM1 

GROW7 


BANDM 


I  DfW.1 


STANDATA 


SCSDI 

GFRDPNC 

BEACONFD 


STAIR 

SCRS8 

SEBA 


SHELL 

PILOT4 

SCFXM2 


SCSD6 

GROW15 

SHIP04S 


FFFFF800 

GANGES 

SCFXM3 


SCTAP2 

GROW22 

SHIP04L 


PDjOTWE 

SIERRA 

SHIP08S 


SCTAP3 

SHIP12S 

25FV47 


SCSD8 

NESM 

CZPROB 


PILOTJA 

SHIP08L 

SHIP12L 


80BAU3B 

PILOT 


Phase  1  Its.  Phase  2  Its 


4 

19  I  93 


Table  2.3  Results  with  SCALE  OPTION  2  and  PARTIAL  PRICE  1 


Analyai »  of  LP  Teat  Problems 


tie* 

era 

ca 

I  Its.  1 

1  Phase  2  I 

Total  Its.  I  Blocked  Its. 


Total  Time  I  Solution  Time 


TaMe  2A  Resalts  with  SCALE  OPTION  2  sad  PARTIAL  PRICE  10 


13664.8 

79939.9 


WTO 


©Sow 
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1 Problem  Nam 


ADUTTLE 

SC205 

SCAGR7 


RECIPE 
VTPBASE 
SHARE  IB 


SCORPION 

CAPRI 

SCAGR23 

SCTAP1 

BRANDY 

ISRAEL 

ETAMACRO 

SCFXM1 

GROW7 

BANDM 

E226 

STANDATA 

SCSD1 

GFRDPNC 

BEACONFD 

STAIR 

SCRS8 

SEBA 

SHELL 

PILOT4 

SCFXM2 

SCSD6 

GROW15 

SHIP04S 

FFFFF800 

GANGES 

SCFXM3 

SCTAP2 

GROW22 

SHIP04L 

PUjOTWE 

SIERRA 

SHIP08S 

SCTAP3 

SHIP12S 

2SFV47 

SCSD8 

NESM 

CZPROB 

PILOTJA 

SHIP08L 

SHIP12L 

80BAU3B 

PILOT 


•i 


_ 0 

202 

16S 

72 

74 

265 

_ 66 

279 

95 

306 

77 

510 

590 

186 

0 

15 

906 

552 

857 

312 

0 

15 

587 

677 

17 

432 

42 

2634 

970 

1261 

619 

1907 

12 

44 

1834 

13352 


Total  Time  I  Solntioa  Time 


2.6 

11.6 

22.1 


896.2 

340.3 
5830.0 


1616.5 

1286.4 

932.7 


5961.1 


13325.9 

86119.7 


364.6 

429.4 

772.2 


431.5 

446.4 

61.9 


3291.2 

734.3 

113.6 


814.8 

263.1 

5763.8 


1541.8 

1180.9 
833.3 


246.7 

579.2 


13073.7 

85900.0 


Table  IS  Results  with  SCALE  OPTION  2  and  PARTIAL  PRICE  20 
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blem  Na 

C5 

1  Phase  1  Its. 

IPhase  2  I 

iBsancai 


Solatioa  Time 


FFFFF800 

GANGES 

SCFXM3 


SCTAP2 


I  Tr;  ;(« 


SHIP04L 


PILOTWE 


SHIP08S 


SCTAP3 

SHIP12S 

25FV47 


SCSD8 

NESM 

CZPROB 


PILOTJA 

SHIP08L 

SHIP12L 


80BAU3B 

PILOT 


0 

33 


492 

610 

17 


430 

43 

2162 


929 

1280 

744 


2032 

12 

46 


1303 

11660 


12000.6 

71991.3 


Table  24  Resalts  with  SCALE  OPTION  2  aad  PARTIAL  PRICE  40 


11747.8 

71772.9 
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Scaled  PP  1 
Objective  Valae 

-4.64753 142857 14E+02 

ADUTTLE 

SC205 

2.2549496316238E+05 
•5.220206121 1707E+01 

SCAGR7 

-2.3313898243310E+06 

-4.1573224074 142E+02 

RECIPE 

-2.6661 600000000E +02 

VTPBASE 

SHARE1B 

BORE3D 

1.2983 146246136E +05 
-7.658931 8579 186E+04 
1.373080394208 5E+03 

SCORPION 

CAPRI 

SCAGR23 

1.878 124822738 1E+03 
2.6900129137682E+03 
-1.4753433060769E+07 

SCTAP1 

BRANDY 

ISRAEL 

1.4122500000000E+03 
U185098964881E+03 
-8.9664482 186305E+05 

ETAMACRO 

SCFXM1 

GROW7 

-7.5571521856065E+O2 
18416759028349E+04 
+1.778781 1814712E+07 

BANDM 

-1.586280184 5012E +02 

-1.875I929066371E+01 

STANDATA 

1 257699 5000000E +03 

SCSDI 

GFRDPNC 

BEACONFD 

8.6666666743334E+00 

6.902235999S488E+06 

3.3592485807200E+04 

STAIR 

SCRS8 

SEBA 

-2.51266951 19296E+02 
9.0429998618888E+02 
1.571 1600000000E +04 

SHELL 

PIL0T4 

SCPXM2 

1 208825 3460000E +09 
-2.5811397876813E+03 
3.6660261 564999E+04 

SCSD6 

GROW15 

SHIP04S 

5.050000007 8262E +01 
- 1.06870941 29358E+08 
1.79871 470044 54E +06 

FFFFF800 

GANGES 

SCFXM3 

5.5567959 102690E+05 
-1.0958576909180E+05 
5.4901254549751E+04 

SCTAP2 

GROW22 

SHIP04L 

1.7248071428571E+03 
-1. 608343364825 6E +08 
1.7933245379704E+06 

PILOTWE 

SIERRA 

SHIP08S 

-2.7200979898204E+06 
1 .5394362 1 83632E+07 
1 .9200982 105346E+06 

SCTAP3 

SHIP12S 

25FV47 

1.4240000000000E+03 
1.4892361 344061E+06 
5.501 8458882867E+03 

SCSD8 

9.049999999254 6E +02 

1.407 6086885553E+07 
2.1851966988566E+06 

PILOTJA 

SHIP08L 

SHTP12L 

-6. 1 1 3062552004 6E +03 
1.90905521 13891E+06 
1 .470 1 879 1 93293E+06 

80BAU3B 

PILOT 

9.8722925484290E+05 

-5.5740380063326E+02 

Scaled  Scaled  Scaled  Uascaled  Descaled 
PP  10  PP  20  PP  40  PP  1  PP  10 


S.9E-10  I  1.1E-09  I  6.4E-10  I  7.1E-09  I  7.5E-09 


3.5E-13 


3.4E-06  3.4E-06  3.4E-06 


1.7E-09  2.9E-10  4.8E-05  3.9E-05 


4.0E-08 

3.4E-08 


3.6E-08 

6.2E-08 


6.7E-07 

4.6E-05 

3.0E-07 

3.0E-07 

1.9E-06  2.4E-07 


2.1E-07 


E-04 

52E 

•04 

E-04 

-06 

6.1E-07 


9.4E-08 


6.1E-07 

4.7E-06 


2.3E-06 


4.9E-07 

1.1E-06 

3.0E-10 

4.7E-04 

1.7E-07 

5.1E-06 

6.5E-05 

Table  3.1  Relative  Error  la  Objective  Compared  to  SCALE  OPTION  2  aad  PARTIAL  PRICE  1 


TO 


•SWi 


TO 


Figure  3.2 

Solution  Time  Factors  Compared  to  No  Scaling,  Full  Pricing 


Figure  3.3 

Solution  Time  Factors  Compared  to  No  Scaling  Equivalents 
Sorted  by  Factor  for  Full  Pricing 


Figure  3.4 

Solution  Time  Factors  Compared  to  Scaling,  Full  Pricing 
Ordered  by  Factor  for  Partial  Price  10 


Figure  3.5 

Ratio  of  Iteration  Count  to  Number  of  Rows 
Ordered  by  Number  of  Rows 
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Picture  for  Problem  1  AFIRO  Magnified  10.0 
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Picture  for  Problem  2  ADLITTLE  Magnified  4.0 
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Picture  for  Problem  3  SC205  Magnified  2.0 
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Picture  for  Problem  5  SHARE2B  Magnified  4.0 
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Picture  for  Problem  6  RECIPE  Magnified  2.0 
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Picture  for  Problem  7  VTPBASE  Magnified  2.0 
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Picture  for  Problem  10  SCORPION  Magnified  1.0 
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Picture  for  Problem  16  ETAMACRO  Magnified  0.5 


Picture  for  Problem  17  SCFXM1  Magnified  1.0 
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Picture  for  Problem  20  E226  Magnified  1.5 


Picture  for  Problem  21  STAND ATA  Magnified  0.5 
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Picture  for  Problem  25  STAIR  Magnified  1.0 


Picture  for  Problem  26  SCRS8  Magnified  0.25 


Picture  for  Problem  27  SEBA  Magnified  0.5 


Picture  for  Problem  28  SHELL  Magnified  0.25 
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Picture  for  Problem  29  PIL0T4  Magnified  0.5 


Picture  for  Problem  30  SCFXM2  Magnified  0.5 


Picture  for  Problem  31  SCSD6  Magnified  0.25 


I 


Picture  for  Problem  33  SHIP04S  Magnified  0.25 
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Picture  for  Problem  34  FFFFF800  Magnified  0.5 
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Picture  for  Problem  35  GANGES  Magnified  0.25 
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Picture  for  Problem  38  GROW22  Magnified  0.5 


AruUfffu  of  LP 


Picture  for  Problem  42  SHIP08S  Magnified  0.125 
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Picture  for  Problem  43  SCTAP3  Magnified  0.125 


Picture  for  Problem  44  SHIP12S  Magnified  0.125 


Picture  for  Problem  45  25FV47  Magnified  0.25 


Picture  for  Problem  46  SCSD8  Magnified  0.125 
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Picture  for  Problem  51  SHIP12L  Magnified  0.0625 


Picture  for  Problem  52  80BAU3B  Magnified  0.05 


Picture  for  Problem  53  PILOT  Magnified  0.125 
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Abstract 


A  set  of  linear  programing  test  problem  is  analyzed  with  MINOS, 
Version  5.1.  The  problem  have  been  run  with  different  options  for  scaling 
and  partial  pricing  to  illustrate  the  effects  of  these  options  on  the 
perforaance  of  the  simplex  method.  The  results  indicate  that  the  different 
options  can  significantly  improve  or  degrade  the  performance  of  the  simplex 
method,  and  that  these  options  must  be  chosen  wisely. 

For  each  problem,  a  picture  of  the  nonzero  structure  of  the  matrix  A 
is  also  presented  so  that  the  problem  can  be  classified  according  to 
structure. 
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